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ABSTRACT 



The Backlund transformation for the three-level coupled Schrodinger-Maxwell equation 
is presented in the matrix potential formalism. By applying the Backlund transformation 
to a constant electric field background, we obtain a general solution for matched pulses (a 
pair of solitary waves) which can emit or absorb a light velocity solitary pulse but otherwise 
propagate with their shapes invariant. In the special case, this solution describes a steady 
state pulse without emission or absorption, and becomes the matched pulse solution 
recently obtained by Hioe and Grobe. A nonlinear superposition rule is derived from the 
Backlund transformation and used for the explicit construction of two solitons as well as 
nonabelian breathers. Various new features of these solutions are addressed. In particular, 
we analyze in detail the scattering of "invertons", a specific pair of different wavelength 
solitons one of which moving with the velocity of light. Unlike the usual case of soliton 
scattering, the broader inverton changes its sign through the scattering. Surprisingly, the 
light velocity inverton receives time advance through the scattering thereby moving faster 
than light, which however does not violate causality. 

^ E-mail address; qpark@nms. kyunghee. ac.kr 
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1 Introduction 



The nonlinear interaction between radiation and a multi-level optical medium has received 
considerable interests for many years. Recently, this topic has attracted more attention 
in the context of lasing without inversion |l|, |^ and electromagnetically induced trans- 
parency (EIT) 0. EIT is a technique for rendering an otherwise optically thick medium 
transparent to a weak probe laser by coupling the upper level coherently to a third level 
by a strong laser field. The transparency for pulses propagating through an optically 
thick medium has been known earlier, particularly for a soliton (27r pulse) in the shape of 
the hyperbolic-secant type and through a three-level medium More recently, there 
appeared exact analytic solutions for a pair of solitary waves, so called matched-solitary- 
wave pairs (MSP), propagating through a three-level medium whose invariant shapes are 
more general than the hyperbolic-secant type [^]. 

In this paper, we present a new type of analysis for the three-level coupled Schodinger- 
Maxwell equation based on the matrix potential formalism. In particular, we find the 
Backlund transformation of the coupled Schodinger-Maxwell equation in terms of the 
matrix potential variable and apply the Backlund transformation to the constant electric 
field background to obtain new MSP solutions which generalize the result in Ref. 0. 
These solutions in general describe the breakup of a MSP into another MSP with slower 
velocity and a soliton pulse moving with light velocity, or the reverse process of fusing a 
MSP and a soliton into another MSP. With a specific choice of parameters, these solutions 
reduce to the MSP solution in in Ref. P] which describes steady-state propogation of MSP 
without breakup. The generality in the shape of a MSP solution is explained through the 
SU (2) X U{1) group symmetry of the three-level A-system with equal oscillator strengths. 
We show that the Backlund transformation also allows a nonlinear superposition rule for 
solitons as well as MSP solutions. An explicit formula for the nonlinear superposition 
is given in terms of matrix potentials and used to generate two-soliton and nonabelian 
breather solutions. We consider in detail the scattering of a specific type of solitons which 
we call "invertons". Invertons are a pair of different wavelength solitons one of which 
moving with the velocity of light and the other with slower velocity. Unlike the usual 
case of soliton scattering, the broader inverton changes its sign during the scattering 
process. Surprisingly, the light velocity inverton receives time advance through the scat- 
tering thereby moving faster than light. We show that however causality is not violated. 
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A typical nonabelian breather describes a breathing Ovr-pulse Ei {v < c) which afterwards 
transfers to the nonbreathing E2 pulse moving with light velocity. 



2 Matrix potential formalism of the A-system 



Consider a A-configuration where level three is higher than level one and two. The system 
of equations governing the propagation of pulses are given by the Schrodinger equation 



dci 
'dt 

dC2 

'dt 
dt 



1^203 



and the Maxwell equation 



dx c dt 
dx cdt ^ 



llCiC^ 

/2C2C3. 



(2) 



Here, /j = 2TTNfifuJi/h; i = 1,2 and c^; k = 1,2,3 are slowly varying probability ampli- 
tudes for the level occupations, Qi = fiiEi/2h are the Rabi frequencies for the transitions 
i ^ 3, El and E2 are the slowly varying electromagnetic field amplitudes, /ij is the dipole 
matrix element for the relevant transition and Ui is the corresponding laser frequency, and 
is the density of resonant three-level atoms. For brevity, we introduce a new coordinate 
z = t — x/c, z = x/c so that d = d/dz = d/dt, d = d / dz = d / dt + cd/dx. 

As in our earlier papers 0, §, the main tool of our analysis will be using a matrix 
potential g instead of the probability amplitudes q in the following way; let (7 be a 3 x 3 
unitary matrix whose second row is the complex conjugation of probability amplitudes, 
i.e., 
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C'i C'2 C3 I , (3) 

y * * * 

where the first and the third rows are to be determined later. In terms of g, the density 
matrix p whose components are pmk = CmCl takes a simple form. 



r9~'Tg, 



T 



/O 


Vo 



0\ 
-ih 
0/ 



(4) 



The specific choice of the matrix T is not essential. One may consider an arbitrary 
diagonal matrix T to handle density matrices in a more general context. Note that the 
first and the third rows of the matrix g do not affect the density matrix p. In other words, 
the density matrix p is invariant under the left multiplication of g by any matrix h, 

9^9= hg, (5) 
which commutes with T and thus of the form 



h= * \ . (6) 




At first sight, introducing the matrix g with more redundant components than q's may 
seem an unnecessary complication. However, this is not so. In fact, it not only manifests 
the symmetry group structure of the system, but it also simplifies the problem of solving 
differential equations. Later, we show that the Backlund transformation of the system, a 
solution generating technique, also takes a simple form in terms of g. The main advantage 
of using g is that g solves the Schrodinger equation identically. In order to see this, we fix 
the redundancy introduced by Eqs. (|) and (^). Adopting the notation for the following 
matrix decomposition; 



[0 
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(7) 



we fix the redundancy by imposing the constraint condition on g, 

{9-'dg)H = 0, {dgg-')H = 0. (8) 

One can always solve the constraint by finding an h which makes g to satisfy the constraint 
via the transform in Eq. (^.0 We also parameterize the remaining components of g^^dg 
such thatR 

' -zfii 
g-^dg = I -tn^ I . (9) 
-iQl -iQ; 



■^The existence of such an h can be proved by adopting the field theory formulation of 
the problem as in fel. However, we do not need the explicit expression of h. 



*Note that g ^dg is anti-Hermitian and g is unitary. 
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The nonzero components in Eq. express the Rabi frequencies in terms of g. In this 
new parametrization, one can readily check that the Schrodinger equation ([T|) arises from 
the simple identity, 

dg^ = dg-' = -g-'dgg-' = -g-'dggl (10) 

This situation may be compared with the ordinary electromagnetism where the static 
electric field E in terms of a scalar potential (p, E = — V0, solves the curl-free condition 
V X ^ = 0. Likewise, we solve the Schrodinger equation in terms of a matrix potential g 
and express the electric field components Qi in terms of g as in Eq. (P). Originally, the 
Schrodinger and the Maxwell equations are made of five complex, component equations 
in total. Three of them (the Schrodinger part) are now solved identically in terms of 
g, where g is partially constrained by Eq. (^). Then, the remaining degree of g can 
be parametrized by two unknown complex functions and the Maxwell equation changes 
into two component equations for these two unknown functions only.0 In this way, the 
Maxwell equation resembles the Poisson equation in electrostatics. However, we do not 
need explicit component expressions in this paper so we supress them. What we need is 
the expression of the Maxwell equation in terms of g such that 

dig-'dg) = Q[T, g-'fg]Q (11) 

where 



T 



/-f 0' 




-1 I , g = k/h I . (12) 

V Of 



Note that this is indeed consistent with Eq. The Maxwell equation in Eq. ([11 

possesses a symmetry under the change; g g = gh where h is an arbitrary constant 

diagonal matrix. That is, g again satisfies Eqs. (§) and (pUf) . This symmetry becomes 

enhanced to a larger one when the oscillator strengths are equal { li = I2 = s) so that 

Q is the identity matrix. For the A-system, the oscillator strength s is positive which 

^ In mathematical terms, we have associated the density matrix p with the coset G j H — 
{SU{3)/ SU{2) X U{1)) and introduced the matrix g for the parametrization of G/H. The 
constraint equation on g makes a specific choice for each equivalence class, and variables 
which parameterize equivalence classes are determined by the Maxwell equation. Since the 
constraint restricts the subgroup H{= SU{2) x [/(l))-part of the variable g, the remaining 
unconstrained part of g can be expressed in terms of two unknown complex functions ipi 
and ip2- In the gauged sigma model formulation, the Maxwell equation is vector gauge 
invariant so that it decouples from the _ff-degree of freedom, i.e. it reduces to two complex 
equations only in ipi and ip2- 



5 



we assume throughout the paper. With equal oscillator strengths, h can be a constant 
unitary matrix of the form 

//in h^2 \ 
h=\h2i . (13) 

V /i33; 

In other words, /i is a constant matrix belonging to the group SU(2) x U{1). In terms of 
physical variables, this symmetry amounts to the transform. 







h* 

'hi 
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(14) 



where the primed variables are solutions of Eqs. (^ and (|TT|) provided that unprimed 
variables are. In particular, if the unprimed solution is a single 27r-pulse ( fii = 0, ~ 
sech(t — x/v)/tp), the primed solution represents a simulton solution [0. When the 
oscillator strengths are equal {Q = 1), the theory becomes integrable and exact analytic 
solutions can be found. However, even for Q ^ 1, the Maxwell equation admits a Lax 
pair representation and infinite sequences of conserved integrals can be found . In this 



paper, we assume the equal oscillator strengths so that the Maxwell equation becomes 

d{g-'dg) = [T, g-'fg]. (15) 



3 Backlund transformation 

The Maxwell equation in Eq. (^) is equivalent to the consistency condition: [Li, L2] = 
of the overdetermined linear equations, 

Li* = {d + g-^dg + XT)-^ = 

L2^ = {d+jg-'fg)^ = 0, (16) 

where A is a spectral parameter. We may apply the inverse scattering method to Eq. 
(|T^) and obtain exact solutions as in ||^. Instead, we present in this paper an alterna- 
tive, simpler method - the Backlund transformation (BT) - which allows a more direct 
construction of exact solutions. Moreover, the advantage of using the matrix potential 
g becomes clear when the BT is used to establish a nonlinear superposition rule for a 
couple of single pulses. This is similar to the linear case where the application of the 
usual superposition rule is easier in terms of the scalar potential rather than the electric 
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field itself. Let go and \E'o be a particular solution of Eqs. (^) and (|16|), then gi is also a 
solution of Eq. if it satisfies the Backlund transformation of type I (type I-BT): 

Type I - BT ; g^'dgi - g^'dgo - Hgi'go, T] = 

tvdig^'9o)+9i'fg,-g,'fgo = 0. (17) 

Here, rj is an arbitrary parameter of the transformation. Once again, the type I-BT is a 
set of overdetermined first order partial differential equations whose consistency requires 
that go and gi should be both solutions of Eq. ([Tsl) . An equivalent expression of the BT 
is in terms of the linear function \I' as in [|1^]. We define the Backlund transformation of 
type II (type II-BT) by: 

Type II - BT ; = ^-(l + ^g-igo)^^. (18) 
A — IT] A 

It can be readily checked that gi and \E'i satisfy Eqs. (|^) and ([T6|) provided Eq. ([T7| ) 
holds and vice versa. The type II-BT is particularly useful for establishing a nonlinear 
superposition rule. Assume that ((?„, \E'a) and {gt, "^b) are two sets of solutions with BT 
parameters rja and rjb respectively, solving the BT for a particular solution, {go, \I'o). If 
we apply the BT once more to {ga, ^E'a) with t] = rjh and also to {g^, "^t) with r] = rja, and 
require that they result in the same solution (this amounts to the commutability of the 
diagram in Fig. 1), then we obtain from the type II-BT, 

^ ^ ^l + ^^-V)(l + T^?.-%)^o. (19) 



(A - iria) (A - irjb) A A 
Equivalently, we have the nonlinear superposition of two solutions, 

9 = {Vb9a - Va9b)9o\Vb9b^ - Va9a^y^- (20) 

Combining this expression with the type I-BT, we obtain a useful formula for the nonlinear 

superposition of the Rabi frequencies, 

1 i 
9~^dg = -{ga^dga + g^^dgb) + -[g'\Va9b + Vb9a), T], (21) 



Or, 



2 



I b ^ 

^2 = 2'^^2 + ^2) ~ 2"^^^ 



F = iVb9b ^ - Va9a ^)9oiVb9a - Va9b) \Va9b + Vb9a)- (22) 
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Figure 1: Commutability diagram for the nonlinear superposition rule. 

4 Matched pulses 

Now, we construct solutions by integrating the type I-BT directly. We choose the partic- 
ular solution go for a constant electric field. 



9o ^9 go 



/ 









—iVti 




-iVLq 

-iVLo 



A 



(23) 



for a real constant fin, so that 



9o = e 



/ cos(r2o-z) — isin(fio-2) 











y— isin(fio-2) cos(fio2) 



(24) 



If we set / = e '^^(/i, the type I-BT becomes 

^r^Br' + r'ff - f 

Since f~^df is anti-hermitian, Eq. (|25|) requires that 

if-' -f,T] = 0, 

which we solve by taking 

f-^- f = 2i sin e 




0. 



(25) 
(26) 

(27) 
(28) 



for an arbitrary real parameter sin^. If we rewrite / in terms of another matrix P, 

f = e~^'gi=cose{2P-l)-tsme, (29) 



8 



then Eqs. (^) and ( pQ]) imply that P is a hermitian projection operator, i.e. 

In terms of P, Eqs. ( p5|) and (|26|) can be written by 

{1- P){d + A-ifiT)P = 
{I- P){ifld -T)P = 0, 

where fj = e'^^rj. Since P is a projection operator acting on the three-dimensional 
we may express P using a three-dimensional vector s = (si, S2, S3) by 

3 

Pij = SiS*/C^SkSl) 
k=l 

which transforms Eq. ( |3lD into a linear one 

{d + A-ifiT)s = 
{ifjd-T)s = 0. 

Since A and T commute, this may be integrated immediately resulting that 

^ i -_ 

Si = Yl [6xp A]ifcUfc, for A = (ir/T - A)z - -Tz 

where u = {ui,U2,us) is an arbitrary complex constant vector. Explicitly, 

51 = (cosh H ^= sinh \^z)ui + ° (sinh y/Kz)u3 

2v K \'K 

fiz sz 

52 = exp(— -)U2 

2 7] 

53 = -^(sinh \^z)ui + (cosh ^/Kz ^= sinh ^/Kz)us, 



where K = rj'^/A — Qq. Finally, Rabi frequencies are given by 

3 

Qi = Qq — 2iri cos 9sisy(^ Sksl 



k=l 

3 

^2 = -2i?7 COS ^S2S3/(^ SfcS^) 

k=l 
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and probability amplitudes are obtained through Eqs. (|^) and (p9D, 

3 

ci = 2 cos OsisyC^Sksl) 

k=l 

3 

C2 = cos9\2s2sl/ (y^ SkSp — 1] + i sin 6* 
fc=i 

3 

C3 = 2cos6's3s;/(^ (37) 
fc=i 

Note that this solution, combined with the symmetry transformation in Eq. (11^), repre- 
sents a rich family of single pulse solutions. For the vanishing {Qq = 0), it becomes 

Qi = —{—2iri cos 6uiulexp[iriz sin 6]) 

1 s s 
0,2 = —(—2iricos6u2Unexp[i(riz-\ — z)sm6 ^cos6']) (38) 

A* T] 1] 

and 

•.osOuiUnexpUrjz z) cos 6' — i-zsinO]) 

7] T] 
X 5 

C2 = — {—\ui\'^ expWjz COS 6 — iO] + |M2p exp[(?72; — 2—z) cosO + i&\ — \u2,\'^ ex\i\—rjz cos 6 

T] 

5 5 

:os 6''U3'U2 exp [ — ^cos6' — i{'riz H — z) sin 6*]) 
r] 7] 



2s 

N = \ui\'^ explriz cos 9] + \u2\'^ expWjz cos 9 zcos9] + \us\'^ exp[—riz cos 9]. (40) 

7] 

This solution for ui 7^ 0,^2 7^ 0, rj cos 9 > describes the transfer of the 27r-pump pulse 
in the limit x ^ —00: 





1 


Cl = 


N 




1 


C2 = 


N 




1 


C3 = 


N 


where 





0,1 — > 0, ^^2 — > — ^^7 cos6'sechSie*^^ 

gSi+ie _ g-Si-ie 
Cl 0, C2 — , C3 — > cos6'e~*^^sechEi 

El = cos 9(r]t- ^^K + (f)i), E2 = sm9(r]t- ~^^ x + 02) (41) 

to the 27r-Stokes pulse moving with the velocity of light in the limit x — > 00; 

Qi — ir] cos ^sech [cos 6*77 (t )] exp sin 6*77 (t )], ^2 — * 

Cl ^ 0, C2-^-e-'', Cs^O, (42) 
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where the arbitrary constants = 1,2 determine sohton positions in time and space. 
In the case of ^ = 0, this transfer of 27r-pulse has been given in WM. For mi = or 



U2 = 0, the solution remains as the steady state 27r-pulse given in Eqs. (^) or (|4^ ) 
respectively without causing any transfer of pulses. This steady pulse, in connection with 
the symmetry transform in Eq. ([l^), is the simulton solution [|^. The free parameter 6 
measures the amount of self-detuning of a pulse from the carrier frequency. The 27r-Stokes 
pulse in Eq. ( ^21) is the same as the usual 27r-pulse {9 = 0) but with the carrier frequency 
shifted by the amount 6w = rj sinO. On the other hand, Eq. (^) shows that the 27r- 
pump pulse receives a time independent phase factor exp[i{s sin^^/c?7)x] in addition to the 
shift of carrier frequency. We emphasize that this detuning has nothing to do with the 
frequency detuning of electromagnetic fields from the resonance line. In fact, our system 
is on resonance and thus the parameter 6 measures the self-generated detuning of each 
pulses. The effect of 6* to a single 27r-pulse is to broaden the pulse shape maintaining the 
2tt area of the envelope which is adjusted by the shift of the carrier frequency. Recall that 
due to the symmetry in Eq. (|T^, a more general expression for a single pulse arises as a 



linear mixture of Qi and ^2 in Eq- (^1) which possess a wide range of free parameters. 

If fio 7^ 0, the solution describes pulses more general than the hyperbolic-secant type. 
For the simplicity of analysis, we assume that 6 = and |?7| > 2 1 fig I- We also restrict 
to the parameters; Ui = ri, M2 = T2, M3 = ir^ for real and rewrite Eq. ( pS] ) for the 
notational convenience as follows; 

si = Aexp{\/Kz) + B exp{—VKz) 

2 7] 

S3 = iCexp{VKz) +iDexp{-\^z) (43) 



where r2 is chosen to be y/A^+CP by an appropriate choice of the coordinate origin. The 
coefficients are defined by 

A = - IH ^ ri 

2^ 2^/K' 2^fK 



B = i(l-^K + ^r3 
2^ 2^/K' 2^/K 



6 = — =ri + - 1 i=)r-i 

2^/K 2^ 2^/k' 

n ^0 ^^^^ ^ ^ \ 
D = 7=r-\ -\ — IH i=)r'i 

2VK ' 2^ ^ 2VK ' 
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K = Iv'-^l (44) 

In the limit where x — oo, the solution takes an asymptotic form: 

rjC 

Qi = — fintanhEr, Qo = =sechS/ (45) 

and 

Ci = , ^ sechSr, Co = — tanhEr, = , =sechSr 

= -\ '-^^{x-vit) (46) 

zrjc 

where the velocity vj is 

= =c, (47) 

2s + - 2r/v^ 

which is less than the light velocity c. Comparison of electric fields in Eq. (^) with 
initial populations in Eq. (^) shows that two laser pulses are arranged in the so-called 
counterintuitive order |^. In the a; oo limit, the asymptotic form of the solution is 

nC 

rii = -JlotanhEi. + nf, = ^ ^ sechSj. (48) 

\/ A? + 



and 



where 



ci = , =sechSjr, C2 = — tanhEiT', cq = , =sechSj? (49) 



2s + (r^2 + 2r^v^) 

[X — vpt) — Ac 



vf 



2t]c 
rf + 2r/v^ 



2s + r/2 + 2r/Vir 



1 52 + ^2 

^° = 2^"a^ 

^ (A2 + C2)e^ocosh(2v^(t-x/c) - Ao) + 

In the far past, this solution represents a matched-solitary-wave pair (MSP) moving with 
velocity vi whose invariant shape is not the hyperbolic-secant type. Eq. (|48| ) shows that 
this MSP in general breaks up into another MSP with slower velocity {vp < vj) and a 
soliton pulse moving with light velocity. Figure 2 shows explicitly this breakup behavour 
with parameters, s = 1, ri = 0.7, = 1, Qq = 0.5 and r] = 1.5. 
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By changing the sign of rj, one could equally consider the reverse process of fusing a 
MSP and a soliton into another MSP. In the specific case where K = tj'^/A — Qq vanishes, 
this solution reduces to the steady state pulse which agrees with the MSP solution by 
Hioe and Grobe 0, 

f2i = — f2otanhS5, ^2 = V^.Qosech'Es (51) 



and 



where 



1 i 
Ci = ^psechS^, C2 = — tanhS^, C3 = ^^sechS^ (52) 
V 2 V 2 



x-vst), vs = 2n'J{s + 2nf^). (53) 



2fioc 

Thus, our solution generalizes the result of Hioe and Grobe and shows that MSP prop- 
agates steadily only when the velocity vs of MSP is specifically given by the parameter 
Qq as in Eq. (|53|). Otherwise, it breaks up as explained before. In other words, in the 
case of MSP where K 0, the hyperbolic-tangent pulse Qi induces a partial transfer 
of the other soliton pulse ^2 into the light velocity one. This is similar to the behavior 
of adiabaton moving with a slower velocity which also generates a light velocity solitary 
pulse during its formation [jlll.0 where the most general expression of our MSP solution 
arises when we incorporate the SU (2) x U (1) symmetry transform as in Eq. (0) as well as 
the self-detuning effect (^ 7^ 0). This solution contains a large number of free parameters 
which account for the variety of pulse shapes and initial conditions required for the atoms. 
The stability of this MSP solution against small perturbations is an important issue W5\ 
which we plan to consider elsewhere. 



5 Invertons and breathers 

The nonlinear superposition rule given in Eqs. (|20D-(p2D allows the superposition of two 
MSP solutions Qa and gb whenever they are obtained from go with BT parameters T]a and 
rjf, respectively. In general, ga and gb could possess different set of free parameters and the 
superposition rule requires only that they are obtained from the same go through BT. If 
ga and gb have same set of parameters except for BT parameters, we could analytically 
^We thank Referee for raising this point. 
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continue the BT parameter in Eq. pO| ) to the complex case in such a way that 

r]a = exp[i9B], Vb = expl-iOs] (54) 

for a certain real parameter 6b, and that g is still unitary and becomes a nonabelian 
breather solution. As an example, we take the one-soliton to be the MSP in the previous 
section with a set of parameters 

Qq = 0, ui = 1, U2 = 1,^3 = —i. (55) 

Then, the nonabelian breather solution has an asymptotic form for x — oo, 

rii = 

sin26'B(exp[(2x/c — t) cosOb] sm^Os + tsinOB) + exp[(t — 2x/c) cosOb] sm{9B — tsmOs)) 



and for x oo, 



cos^ 6b cos(2t sin 6^^) — cosh[(4x/c — 2t) cos 9b] sin^ 9b 



sin2^B(e^c; sin(0B - Zs) + e"^^ sin(^B + Zs)) 



—1 + cos^ 9b cos{2Zs) — cosh.{2Zc) sin^ 9b 

^2 = 

Zc = {t — x/c) cos 9b, Zs = (t — x/c)sm9B- (57) 

This describes a breathing Ovr-pulse ^2 {v < c) which transfers to the nonbreathing Qi 
pulse moving with the velocity of light. This is a typical motion of a nonabelian breather. 
Figure 3 shows a breathing motion with 9b = 1-2. Different values of parameters Ui in 
general distort the shape of Ovr-pulse so that the time area of the pulse is nonzero, but 
they lead to the essentially same type of transferring motion except for the case ui = 
or ^2 = where it becomes a steady state breather. 

For real BT parameters r]a and r]b, the superposed solution g in general describes the 
scattering of two MSP solutions. Here, we concentrate only on a particular case of two 
soliton solutions, which we name as a pair of "invertons" due to their interesting new 
features. Consider one solitons ga and gb given by 

/-I 

ga = cos 0a siu 0„ 

\ sin (pa — cos 0, 
COS 0a = tanh[?7at x] 



cr]a 

Tj^ -\- S 

sin 0a = sech[77at x] Vb't 

CVa 
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(56) 



and 

/ cos0fc sin 05 
Qb ^ 0-1 

\ sin 05 — cos 05 
cos 05 = tanh[?75(t — x/c)] 

sin 05 = sech[7;5(t — a;/c)]. (59) 
In terms of Rabi frequencies, these correspond to 

= 0, Q^"^ = -i77„sech[77„t - ^^^^x] (60) 

and 

nf) = -i?75sech[?76(t - x/c)], Q^''^ = 0, (61) 

that is, they represent two 27r-pulses with different resonance frequencies, one moving 
with hght velocity and the other moving with slower velocity. 

The nonlinear superposition rule gives rise to the superposed Rabi frequencies and 
probabilities such that 

ril cos 0a + r]ar]h{l + cos 0a) + 775 



Qi = — i?75sin( 



'^a + ^6 + ^a^6(l + COS 0a + COS 05 - COS 0a COS 05) 



O A Va + V^Vbjl + COS 05 + ?75 COS 05 , . 

ih = -Z77aSm0a^- — — — --r (62) 

'nt + 11b+ VaVbi^ + cos 0a + COS 06 - COS 0a COS 06) 



and 



Cl = 

C2 = 



{Va + Vb)VbSin 0a sin 05 



Vl + Vb+ VaVbi'^ + COS 0a + COS 06 - COS 0a COS 06) 
-jvl + Vb) COS 0a - rjaribi^ + COS 0a - COS 06 + COS 0a COS 06) 
+ ^6 + ?7a?76(l + COS 0a + COS 05 - COS 0a COS 06) 



^ -(^g + VaVbi'^ + COS 06) + ryg cos 05) sin 0a 

^ + ^6 + ^a?76(l + COS 0a + COS 06 - COS 0a COS 06) ' 

This solution describes the scattering of two invertons where the fast moving inverton Ip 
passes through the slow moving one Is- Before the collision, their asymptotic forms are 
given by the following configurations; 

Inverton Ip ; cos 0a = 1 

Cl = 0, C2 = -1, C3 = 

fli = —ir]bsech.[r]b{t — x/c)], 1^2 = (64) 
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and 



Inverton Is ; cos (j)^,— ! 



ci = 0, C2 = — tanh[77ai x], C3 = — sechfT^at ^ x] 

cria cr]a 

Vt^ -\- s 

Vti = 0, 0.2 = -irjasechhjat - — x]. (65) 

cria 



After the collision, their asymptotic forms are 
Inverton Ip ; cos 0a = —1 



Ci = 0, C2 = 1, C3 = 

^1 — •?^bsign( ^" ^^ )sech[77b(t — x/c) + ^], = (66) 



and 

Inverton /g ; cos 0^ = — 1 



tF' -\- s 

Ci = 0, C2 = — tanh[77ai x + 5], 



cr)a 

Va-Vb' " ' Crja 



C3 = — sign( jsech[ryai x + o\ 



2 

n^ = 0, ^2 = -iriasigni "^" ^ )sech[r)at - ^^^^x + S], (67) 

Va - Vb crja 

where sign denotes the sign function and the phase shift parameter 5 is 

S^ln\!I^±^\. (68) 

Va - Vb 

The appearance of the sign function is a novel feature of soliton scatterings in a three- 
level A system which does not arise in the two-level atom case. A careful analysis of 
the above asymptotic solutions shows that the broader pulse changes its sign after the 
collision (thereby the name inverton) while the sharper pulse remains the same. Also, the 
phase shift always arises in such a way that the slow inverton receives time retardation 
while the fast one recieves time advance. This implies that the fast inverton moving with 
hght velocity moves faster than light through the scattering process! However, this does 
not violate causality. In fact, the fast inverton has infinitely stretched tails which trigger 
the slow inverton to transfer its energy to the fast one. 

Figure 4 shows the scattering of the invertons with parameters rja = 0.5, rjb = 0.9 and 
s = 1. In order to see that causality is not violated, we consider the energy conservation 
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law given by 

d d 

where 

p = + 1^2!' + |ci|' + |C2|' = + l^^aT + l- Icsl' 

J = c{\Q,\' + \Q2\'). (70) 

Figure 5 shows the energy transfer process between the two invertons. When the right 
end tail of the fast inverton Ip reaches the slow inverton Is, the energy is transferred from 
Is to If and Ip receives phase advance. After the center of Ip reaches Is, the energy 
of Ip is absorbed into the Is which results in a retarded phase shift of Is- Thus, no 
energy is really transferred faster than light. The scattering process is indeed a nicely 
balanced exchange process where the scattered invertons maintain their shapes invariant 
except for the change of sign. The change of polarizations of simultons in the scattering 
process has been known earlier |jl6| and the scattering of two single-frequency 27r-pulses 



with different resonance frequencies is also considered in |T^. However, the binary, sign 
changing behavior of invertons is a new aspect which did not receive attention in earlier 
works, and hopefully will find some applications in the future. 



6 Discussion 

In this paper, we have introduced a new formalism for the three-level coupled Schodinger- 
Maxwell equation based on a matrix potential variable, and through the Backlund trans- 
formation we obtained various new solutions which generalize previously known solutions. 
We found new MSP solutions, more general than the steady state MSP found by Hioe 
and Grobe, which account for the generic breakup behavior of a MSP into another MSP 
and a light velocity soliton. Two-soliton solutions and nonabelian breathers are obtained 
through the nonlinear superposition rule and scattering of invertons are analyzed in de- 
tail. In this new formalism, we have demonstrated the SU{2) x U{1) group symmetry 
of the three-level system with equal oscillator strengths and using that we have found a 
general expression of solutions which mixes pulses with different velocities. This group 
symmetry also accounts for the variety of steady state MSP solution by Hioe and Grobe 
which linearly superpose different Jacobi elliptic functions. Another new feature of our 
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formalism is the introduction of a self-detuning parameter 6 to the solutions. For each 
27r-pulses, these ^-values are conserved in time and survives from the scattering process. 
That is, it is a conserved charge which can be used to label each pulses in addition to the 
area of pulses. We also emphasize that our new formalism is not restricted to the three- 
level case but also can be extended to other multi-level cases ^ with the same analytic 
power. 

Finally, we comment on matched pulse propagation through absorbing media. In this 
case, the time rate equation for the probability amplitude C3 in Eq. (0) is replaced by 
dc^/dt — > dcs/dt + 7C3, which implies the decaying of the excited state |3 > at a rate 7 
to states other than |1 > and |2 >. This obviously breaks the probability conservation 
law; |cip + |c2p + |c3p = 1 so that our matrix potential formalism does not apply to 
this case. Nevertheless, it is important to note that the SU{2) x U{1) group symmetry 
of the three-level system with equal oscillator strengths persists even in the case of an 
absorbing medium. This shows that recent analytic and numerical studies on matched 
pulse propagation through absorbing media ||T^-|^ can be extended to more general 
initial conditions. For instance, most studies assume the standard initial condition of the 
population being in the ground state. The SU{2) x f/(l) group symmetry as indicated in 
Eq. (|l^) maps the standard initial condition to a coherently prepared one which linearly 
superposes |1 > and |2 >. This also maps two Rabi frequencies Qi and ^2 via the 
linear transform as in Eq. (|I4|) so that one obtains a more general description of pulse 
propagation. Thus, for example, when a matched pulse of equal amplitude is found, the 
above symmetry generates a set of new solutions with unequal amplitudes. This agrees 
with results in the recent numerical study which showed the existence of stable, shape 



invariant matched pulses of unequal amplitudes depending on the coherent preparation 
of initial conditions. 
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